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Abstract 

Exact integral expressions of the skew orthogonal polynomials involved in Orthogonal 
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rank) skew orthogonal polynomials are average characteristic polynomials of random 
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1 Introduction 



Families of Orthogonal (or skew-orthogonal) Polynomials, have many applications to 
mathematics and physics 

Here, we will have in mind applications to Random Matrix Theory (RMT) 
@)@)H)§|)0]) i- e - disordered solid state physics ||, QCD J7|, or statistical physics on a 
random fluctuating lattice |||8j (2D quantum gravity). In all these fields of physics, one 
is interested in the spectrum of a matrix (Hamiltonian, Transmission matrix, S-matrix, 
Dirac operator,...), which can be considered as random for various reasons (disorder, 
random impurities, quantum fluctuations, chaos or non-integrability,...). It was ob- 
served that the spectrum of a large random matrix shows universal properties |10[ 



(2-point correlation function, in the short or long range regime; universal conductance 
fluctuations of mesoscopic conductors). One possible way to understand and prove 
that universality is through the "orthogonal polynomials" method, which we shall re- 
call below. In order to extract some useful numerical results, it is important to have 
some asymptotics of the orthogonal polynomials in some special limit. 

The type of orthogonal polynomials involved, depends on the symmetry of the 
matrix ensemble 0. The case of a physical system with broken time-reversibility (for 
instance a mesoscopic conductor in the presence of a magnetic field), represented by a 
U(N) invariant matrix ensemble, was extensively studied, because it is the simplest ||. 

Here, we shall focus on the O(N) and Sp(2N) invariant matrix ensembles, which 
appear for physical systems with time-reversibility and/or half-integer spin with broken 
rotational symmetry. These ensembles involve families of skew-orthogonal polynomials. 

The aim of this article is to present a remarkable exact expression of the skew 
orthogonal polynomial as an integral, and deduce from it the required asymptotics. 

Section 2 is a brief introduction to the orthogonal polynomial's method in RMT, in 
section 3 we give and prove the remarkable exact expressions for the skew-orthogonal 
polynomials, and in section 4, we consider their asymptotics. 

2 The Orthogonal Polynomials 

Consider the partition function of a random matrix M: 

zf[V] = [ dMe"^ try(M) (2.1) 

where E$ is the set of all iV x iV real symmetric matrices, E$ is the set of all 
N x N hermitian matrices, E$ is the set of all 2N x 2N self-adjoint real quaternionic 
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matrices^, and dM is the Haar measure on V(x) is a polynomial potential, 

bounded from below, and Np = N, N, N/2 respectively for (3 = 1, 2, 4. 

The angular degrees of freedom of M can be integrated out, and 2A_ can be rewritten 
as an integral over the N eigenvalues (Ai, . . . , Ajv) of M only ]TT|,^]: 



A' 



zf[V] = uf \\d\ % |A(A)[' 



(2.2) 



where U^' is the volume of the group O(N), U(N) or Sp(2N) respectively for (3 
1,2,4. dA = d\e~ NV ^ is the measure element, and 



A(A) = J](A, - A,) 



(2.3) 



i<j 



is the Vandermonde determinant, which can be rewritten as: 

/I Ai \\ 
1 A 2 



A (A) = det 



A^ 1 



Vl A 



N 



X 2 



det(Af) = det Pj(\i 



(2.4) 



where Pj(X) = A J + . . . is an arbitrary monic polynomial of degree j. The last equality 
is obtained by linearly mixing columns of the determinant, and the first equality is the 
well known Vandermonde determinant, which can be found in any math textbook 0. 



The computation of integral |2.2| becomes easier with a special choice of the poly- 
nomials Pj(X), chosen orthogonal with respect to an appropriate scalar product ||1 1|| : 



In the unitary case (3 = 2, the scalar product under consideration is: 



<f\g>- 



dxf(x)g(x) 



and the polynomials P n {x) are chosen orthogonal: 



< P I P >: 



(2.5) 



(2.6) 



In the orthogonal case (3 
symmetric: 

< f\g >=- < g\f >-- 



1, the scalar product under consideration is skew- 



dxdy f(x) sgn(x - y) g(y) 



and the polynomials P n (x) are chosen skew-orthogonal: 



< p2n\P2m > — < P2n+l\P r , 



2m+l 



< P2n+l\P2m > z 





h n 5 n 



2 M £ Ejf can be viewed either as a 2N x 2N matrix with complex number entries or a N x N 
matrix with quaternion entries. It has N eigenvalues, each twice degenerated. [|| 



(2.7) 

(2.8) 
(2.9) 
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In the symplectic case j3 — 4, the scalar product under consideration is skew- 
symmetric too: 



< f\9 >=-<9\f>= dx (f(x)g'(x) - f'(x)g(x)) (2.10) 

J — oo 

and the polynomials P n {x) are chosen skew-orthogonal: 

< P2n\P2m > = < P2n+1 \ P2m+1 >= (2.11) 
< P2n+l\P2m >= h n 5 nm (2-12) 



In all three cases, the partition function 22 reduces mainly^ to Z — Yln=i ^n-i> 



where np is called the "Fermi level" by analogy with a system of fermions: 

n F = ^N„ for = (1,2,4) (2.13) 

When iV is large, most of the physical quantities and relevant observables are related to 
properties of h n in the vicinity of the Fermi level: n — > oo, iV — > oo and n — riF ~ 0(1). 

2.1 Determination of the orthogonal polynomials 

For a generic potential V(x), those orthogonal polynomials exist, and can be con- 
structed by recurrence. Indeed, we start from Pq(x) = 1, then the coefficients of Pi are 
determined by the orthogonality conditions, and by recurrence, we determine P n and 
h n for all n. 

Note that for the skew-orthogonal polynomials, there is an ambiguity: P2 n +i is 
defined only up to an arbitrary linear combination with Pm- If one wants a unique 
definition, an extra condition should be added, for instance that the term of degree 2n 
in P271+1 vanishes. Anyway, the values of h n don't depend on this ambiguity. 

The determination of the orthogonal polynomials by recurrence is unefficient if 
one wants to compute P n for n large. The aim of this article is to present a closed 
expression of P n for any n, and to derive from it some asymptotics in the large n limit, 
and particularly near the Fermi level n — rip ~ 0(1). 

3 An exact expression of the skew-orthogonal poly- 
nomials 

• In the unitary case /3 — 2, it is known that 



3 The actual result may depend on the parity of AT. Details can be found in 



3 



Pi 2) {x) = \ [ dM det(x-M) e ~ NtlV M (3.14) 



f/ (2) 

<-^ra 

Z/ri 



The n th orthogonal polynomial is the average of the characteristic polynomial of 
a n x n hermitian matrix with respect to the weight e - NTrV ( M ) : 

Pi 2 )(x) = (det(x-M)) nxn (3.15) 



This has been known for more than a century |12| (in the context of RMT, see 
e.g. |l|,p|,^3[|). We are now going to generalize this expression to (3 = 1 and 4. 

Orthogonal case (3 = 1. We will prove below that: 



p£\ x ) = _L_ / dM det(x-M) e^^W (3.16) 

^2n ^ A/e - B 2n 



^2n 



y dA 1 ...dA 2n niA i -A,i n^-^ 



7 ( 

^2(1 " i<j 



and 



(det (x-M)) 2nx2n 

^2n+i(^) = -4ty / dM (x + tr M + c n ) det (x - M) e^ 11 y(M) (3.17) 
J dAi . . . dA 2 n |A, - Aj| (2 + ^ Ai + c n ) JJ(x - A*) 



7 (1) 

^2n " i<j 

((x + Tr M + c n ) det (x - M)) 



2nx2n 



the constants c n can be chosen arbitrarily, the choice c n = is the one such that 
the term of degree 2n in P 2 n+i vanishes. 

Symplectic case (3 = 4 



p^(x) = -L / dM det(x-M) e"T try(M) (3.18) 



r/ 1 ' 



n 

(4) /■ _ 

y dAL-.d^JIlAi-A/JJCx-A^ 



7 (4) 



i<3 

(det (x - M)) nxn 



4 



and 

P^ +1 (x) = I dM(x+ trM + c n )det(x- M) e~% tvV{M) (3.19) 

%n J i<j i i 

= ({x+ TrM + c n )det(x-M)) nxn 
3.1 Proof of 

Note that it is sufficient to prove that 

< P 2n \x m >= and < P 2n+1 \x m >= for all m < 2n - 1 (3.20) 
Consider: 

<P 2n \x m > oc J dx dy dAi . . . dA 2n (3.21) 

iica, - \j) H(x - \i) n sgn ^ - x^ sgn ^ - V ) ^ 

then write x = X^n+i- 

<P 2n \x m > oc/ dj/dAi-.-dAan+i y m n (Ai — Aj) (3.22) 

l<i<j<2n+l 
2n 

] [ sgn (Aj - Xj) ] J sgn (Ai - A 2 „+i) sgn (A 2 „+i - y) 

l<i<j<2n+l i=l 

symmetrize with respect to the first 2n + 1 variables: 

2n+l „ 

<P 2n |x m > oc ^ / dy dAx . . . dA 2n+1 y m J] (Ai — Xj) (3.23) 

fc=l ^ l<i<j<2n+l 

] | sgn (A, - Xj) Yl sgn (Ai - X k ) sgn (A fc - y) 

l<i<j<2n+l l<i^k<2n+l 

Note the following identity: 

2n+l 2n+l 2n+l 

IJ sgn (y -Xi) = ^2 s S n (y ~ X k) JJ sgn (A* - Ai) (3.24) 

i=l k=l i=l,i^k 

which gives (and note y = A 2n+2 ): 

< P2n\x m > oc / dAi . . . dA 2n+2 A^ +2 Yli<i<j<2n+i(^ ~ X j) 

rTl<i<i<2n+2 s S n ( X i - X j) ( 3 - 25 ) 



The second bracket is completely antisymmetric in the 2n + 2 variables, so that we 
have to antisymmetrize the first bracket as well. The result is zero when m < 2n, 
because any non-zero antisymmetric polynomial of 2n + 2 variables must have degree 
at least 2n + 1, while the first bracket is a polynomial of degree at most 2n in any of 
its variables. 

By the same argument, one would find that 

< iWi|* m > oc J dAi . . . dA 2n+2 [A™ +2 (c n + E£f A,) ni< i<i < 2ri+1 (A i - A,)' 

rTi<i<i<2n+2 s S n (Ai - Aj)J (3.26) 
which, by antisymmetrization of the first bracket, vanishes when m < 2n — 1. 

3.2 Proof of ffTT3| 

Again, it is sufficient to prove that 

< P 2n \x m >= and < P 2 n+i\x m >= for all m < 2n - 1 (3.27) 



< P2n\x m > oc / da; dAi . . . dA 



(3.28) 



n (A.-A^n^-A.) 

l<i<j<n 



2 ' mi" 14 - x m 



X - A, ; 



Introduce n extra variables (/ii, . . . , /i n ), and consider the 2n x 2n Vandermonde deter- 
minant of the 2n variables (Aj, fii), divide it by n^Aj ~ aO an d take the limit ^ — > A^. 
You get: 



(A, - A,) = km ff?r — - 



With the same trick, we have: 



det 



/l \i A 2 

1 A n A 2 

1 2Ai 

\0 1 2A n 



A 2 "" 1 \ 



A 



2n-l 
n 



2n-2 



(2n- 1)A 2 



(2n-l)Af- 2 / 
(3.29) 



II (Ai-A^JI^-A, 

l<i<j<n i=l 



A ( x, Xi,fii) 



(3.30) 
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fix x 2 
1 Ax X\ 



•'in. 



det 



and 



1 A n 

1 

\o i 



/l Ai A? 



\2 

2Ai 



2A n 



Af 



\2n 

. . 2n\j n - 1 



2nX 



2n-l 



d_ 

dx 



Tl (A,-A J ) 4 n^-A,) 2 = det 



l<i<j<n 



i=l 



2n 



1 




Vo 



A» 
1 
1 



A 



2 
n 

2x 
2Ai 



A 



2» 



(2n + 
2nA? 



2n-l 



2n-l 



2A T 



2nX 



2n-l 
n 



Therefore, the integrand in |3.28| is a (2n + 2) x (2n + 2) determinant: 



det 



x 



fl x 

1 A x X 2 

1 A n X 2 n 

1 2x 

1 2Ai 



x 



2n 



X'' 



\ 



A? 



2d 



A 



2d 



(2n + l)x 
2nA? 



2n-l 



mi 





m— 1 



2n-l 







2A, 



2nA 



2n-l 



o / 



we note x = X n +i, and by antisymmetrization, it becomes: 



det 



fl Ai 

1 A n+ i 

1 

Vo i 



A? 



X 2 

A n+1 

2Ai 



Af 



2d 



K \ 



\2n \m 
A n+1 A n+1 

2nX 2 1 n - 1 mAf- 1 



2A 



n+l 



2n\ 2 n n ~ 1 1 m\ m - 1 



v n+l 



(3.31) 



(3.32) 



(3.33) 



which obviously vanishes when m <2n. 

By the same argument, one would find that < P2 n +i\x m > reduces to the same kind 
of integral, but with m replaced by m + 1, and vanishes when m < 2n — 1. 

We have thus proven that the skew-orthogonal polynomials are indeed given by 



|3~T6|andEn^ 



7 



4 Large N asymptotics 



Most of the large N universal statistical properties of a random N x N matrix M 

pi/3] 



belonging to one of the three ensembles can be expressed in terms of a few h n , 



with n close to the Fermi level 

n F = |iV> (4.34) 
More precisely, for /3 — 2, we need asymptotics of P n in the limit 

iV^oo , n^oo , n-N~0{l) (4.35) 

for /3 — 1, we need asymptotics of Pm and P2n+i in the limit 

iV^oo , n^oo , 2n-N~0(l) (4.36) 

and for (5 = 4, we need asymptotics of P2n and P2n+i in the limit 

iV^oo , n^oo , n- N = n- 2Ni~ 0{l) (4.37) 

4.1 The resolvant 

We introduce the function W(z) usually called the resolvant or Green function: 
where M e E^m and the mean value is taken with respect to the weight: 

e -m,TrV(M) ^gg) 

When there is no ambiguity, we will drop the (3, m or V indices, and write the resolvant 
as W(z). Note that we have chosen a normalization such that: 

W(z) ~ - (4.40) 

The reason to introduce the resolvant is that the logarithmic derivative of P n (x) is 
proportional to the resolvant W m (z) (from |3~T4] , [3 . 1 6j , j3 . 1 8| , at least when n is even) for 



some appropriate value of m, and with a potential of the form: 



V(z) = ^V(z)-rhi(x-z). (4.41) 



More precisely, we have: 
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In the Unitary case (3 — 2: 

^i = nW n (z)\ z=x with V(z) = -^)--ln(x-z) (4.42) 

Pn (X) n n 

i.e. m = n, r = - and T = jj (— > 1 when n — > rip). 
In the Orthogonal case (3 = 1: 



--2niy 2n (2)| z _ a; with V(x) = —V \n(x-z) (4.43) 

P g\ x ) *-* m m 

i.e. m = 2n, r = ^ and T = ^ (— >■ 1 when n — > 

In the Symplectic case /3 = 4: 

2n W - 2n W n (z) \ z=x with V(x) = - V - - ln(x - z) (4.44) 



P 2 ( J(x) ■ ' n n 

i.e. m = n, r = - and T = (— > 1 when n — > rip). 

In all three cases: T = — and r = 

np in 

4.2 Asymptotics for the Resolvant 

In a potential V, the resolvant W(z) = W m (z) satisfies the equations of motion (result- 
ing from invariance of an integral like eq. |2.1| under a change of variable M — > f(M)): 

W(zf - ^W'(z) = \V(z)W(z) - Q{z) + 0{l/n 2 ) (4.45) 
2n p 

where 77 = (1,0, —1) respectively for (3 = (1,2,4) and Q(z) is a polynomial^ of degree 
deg V — 2, which is not determined by the equations of motions, it has to be determined 
by analytical considerations, for instance the one-cut asumption. 
Here, we will consider a potential V of the form: 

V(z) = —V(z)-r1n(x-z) (4.46) 

and we will be interested in the limit where T — 1 and r are small of order 1/n. 
The method is to find first the solution W(z) at T = 1 and r = 0. We write it: 

W(z) = W (z) + ^-Wi(z) + 0(1 /n 2 ) (4.47) 



4 



when V' has poles, Q may have poles too. Q(z) is a rational function, whose poles must be chosen 



in order to cancel the poles of W(z) in eq.4.45 
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and then, add the variations: 

(T- l^Wo+r^Wo (4.48) 

(at order 1/n, we don't need to consider the variations of W\ with respect to T and 
r), the derivatives are taken at T = 1 and r = 0. 

4.3 Contribution of Wo 

The function Wq(z), (as well as its derivatives with respect to T and r) has been 
extensively studied in RMT. Note that Wo is nearly the same for f3 — 1, 2 or 4. Let us 
recall here some of the main features of Wq in order to fix the notations. 



At n — > oo (and T = 1 and r = 0), |4.45| reduces to a quadratic equation for Wo(z). 
The one-cut-solution is: 

Wb(*) = - (V\z) - M{z)^{z-a){z-b)) = iy'(s) - mp(z) (4.49) 

Where M(z) is a polynomial of degree d — 1 (<i = deg V), which is completely deter- 
mined by the large z limit condition [4.40| : 



M(z) = Pol S| (4.50) 
y/(z - a)(z - b) 



The end-points a and 6 too, are determined by [4.40| which implies 



d, = , / , ZV ' {Z) dz = 2W (4.51) 



y/(z -a)(z- b) J y/(z -a)(z- b) 

where the contour encircles the cut [a, b] in the trigonometric direction. 
The imaginary part of the resolvant 

p(z) = ±M(z)y/(z-a)(b-z) (4.52) 

/J7T 

is the average density of eigenvalues of the random matrix (as an example, consider 
the gaussian case: V is quadratic, i.e. V is of degree d — 1, thus M{z) is a constant 
and p{z) = ^ \z — a)(b — z) is the famous semi-circle law). 

Some notations: 

It will be convenient to parametrize z as: 

a + b b - a , b — a ,. Kn . 

z = ^— + cos0 , a := — — (4.53) 

We write: 

ct(^) := (z — a)(z — b) , a/ct(z) = 2ia sin (4.54) 

Note that 0(z) is a multi-valued function. We will see that both determinations 
and —<f) will enter the asymptotic expression of the orthogonal polynomials when 
z E [a, b}. 
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4.4 Variations of Wo with respect to T and r 

It can be proven (see |H| for instance) that 

d<f>(z) 



W T {z) := £=TW (z) = — L=, 

dT V a ( z ) dz 



and 



W r {z) :=P—W (z) 
dr 



1 yJo{z) - \Jo{x)_ 1 



o [Z 



z - X 



In particular at z — x, we have: 



a'(x) 1 



WJx) r - 
2a(x) y/a[x) 

4.5 Contribution of W\ 

For T = 1 and r = 0, and at order 0(l/n), the equation of motion reduces to: 

W 2 (z) - ±W'{z) + 0(l/n 2 ) = jY(z)W{z) - Q{z) 
and we expand W(z) at first order in 1/n as: 



(4.55) 



(4.56) 



(4.57) 



W{z) ~ W (z) + (z) + 0(l/n 2 ) 



At order 1/n, eq. [4.58| gives (using the value of W / o(^) from [4.49| ): 

/3 M(z)v^(i) 



(4.58) 



(4.59) 



(4.60) 



where Qi(z) is a polynomial of degree d — 2. 

Let us factorize M(z) (recall that d = deg V and (7 is the leading coef. of V): 

d-l 

M{z) = g\[{z-z k ) (4.61) 

fe=i 

and decompose W\ in single pole terms. The condition that W\(z) is regular when 
z = Zk allows to determine the polynomial Qi(z), and eventually we get: 



W x {z) 



a'(z 



Aa(z 



d-l 



d 



z - z k )y/a{z) 



fc=l 



(4.62) 



With the parametrization z = 2±2 4- 2a; cos and = + 2a cos 0^, we have 

d-l 

sin 



^ d-l 

-lncr(^) + J^hif 



fc=i 



+ fc d 
"2— } -2 ( 



(4.63) 
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4.6 Asymptotics of the skew-orthogonal polynomials 

Eventually we have computed all the contributions to the asymptotics of the resolvant: 

(4.64) 



W(z) ~ ^W (z) + ^j^Wt(z) + T -W r {z) + ■^W 1 (z) + 0(l/n 2 ) 



i.e.: 



2nW{z) ~ pNpW (z) + (2n - /3Np)W T (z) + W r (z) + rjW^z) + 0(l/n) (4.65) 
where W , W T , W r , W x are given by |44£|, [D>|, |45| (or ^57|) , |46| (or gggj) . 
Combining everything together, we get: 

• P = 2. From P' n jP n = we get the asymptotic orthogonal polynomial (already 



known 14,13, 



^i 2) (z) 



-NiTrf*p{y)dy i(n-N+l)<f, 



+ C.C. 



^/2za sin 

The normalization constant Cn^ = a n+ ^ is such that P n (x) ~ a; n for large x. 



(4.66) 



4.66| is basically the contribution of Wo, which is the same for all three cases /3 — 1, 2, 4. 



The P — 1 and P — A cases contain an extra contributions from Wi. 
• P = l. From P' 2 j p 2n = 2nW(x) we get: 



(i) 



-Niir f* p(y)dy i(2n+l-iV-§) 



v27as!m 



M + (0) 



c.c. 



(4.67) 



where 



d-l 



M+(0) = M_ 



| J 2z sin 



+ 4>k 



(4.68) 



note that M(x) = ga 1 M + ((/>)M_(0), where g is the leading coefficient of V'(x). 



d-l 



-i0 fe /2 



(4.69) 



fe=i 



is the normalization constant chosen so that P2n{x) ~ x 2n for large 

The odd polynomial is found from Pi n ^.\jPi n =< x + Tr M + c n > and 
< Tr M >= 2n\im z __ t00 z 2 {W{z) - 1) (note that we need not |Q7| ). The whole 
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x dependance of Pm+xjPin comes from x + lim 2 ^oo z 2 (W r (z) — 1) = a/ cr(x) — 
Therefore, (and up to an arbitrary linear combination with P2 n ), we have: 



c.c. 



(4.70) 



(3 = 4. From P! l jP 2n = 2nW(x) we get: 




with normalization constant: 



(4.71) 



d-i 



and 




(4.73) 



Note that we have used: ip(x) = -f-a d M + (0) M_(0) 2isin0. 



Some remarks: 

- the derivation presented here is actually valid only when x ^ [a,b], giving only one 
exponential term, with the determination of <f>(x) (from |4.53| ) such that P n (x)e~ NV ^ 
decreases when x — > oo. When x G [a, b], a carefull analysis shows that both deter- 
minations of <p{x) must be taken into account. The only effect is to add the complex 
conjugate exponential (c.c.) to the asymptotics, so that P n (x) is indeed real when 
x G [a, b). Outside [a, b], P n e~ NV decreases exponentially, and in [a, b], it oscillates like 
a cosine function, and it indeed has n zeroes. 

- Our derivation was carried out only in the "one-cut" case. We have assumed that 
the support of p{x) is connected and is made of one interval [a, b\. 

- Those asymptotics are not valid when x is close to a or h. 

- Note that the above expressions all have the correct large x behaviour: P n (x) ~ x n . 
It can be seen easily if one remembers that x ~ ae 1 ^ when i<p — > +oo. 
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4.7 Check of orthogonality 



We have presented a derivation of the asymptotics ( |4.67| - P~T5| ), so that there should be 



no reason to doubt they fulfill the orthogonality condition. However, it is interesting 
to see how. We will just sketch the procedure: 

In all cases, we have to compute integrals of P n P m e~ NV , with x running from — oo to 
+00. The contributions outside [a, b] are exponentially small, the integrals can thus 
be computed inside [a, b]. Within [a, b], terms which oscillate exponentially fast like 
e Nm f p : average to zero at order 0(1/N), so that at leading order, it is sufficient to 
consider only the cross-terms in the product P n P m , with opposite signs for the two 
determinations of <b. 

In the (3 = 1 case, the scalar product < P n \P m > of ^T7\ can be computed by integration 
by part. For that, you need a primitive of P n e~ NV , which is achieved at leading order 
by dividing [4.671 or |4.70| by p(x) = cte M + (0)M_((/>) sin0. 



In the (3 = 4 case, you need a derivative of P n e~^ , which is achieved at leading order 
by multiplying [4.71| or [4.73| by p(x) = cte M + ((f>)M_((f)) sin0. 



Then you find that in both cases ((3 = 1 and 4), and up to unimportant constant 
factors, you have at leading order in 1/n: 

f 11 , sin2(n — m)d) 

< P2n\P2m >0C / d0 \ - ^ = (4.74) 

Jo sin 

/*7T 

< -p2n+i|-P2m+i >oc / d<f> sin(j) sin2(n — m)(j) = (4.75) 
Jo 

/*7T 

< P 2n+ i\P 2m >oc / d0 cos2(n - m)0 oc 5 nm (4.76) 
Jo 

which confirms that our asymptotics indeed fulfill the orthogonality properties. 
Taking into account properly the constant factors, we can determine the h n 's: 

• (3 = 2: 

~ 2vr a 2n+1 (4.77) 

. (3 = 1: 

h$ ~ -^T7 « 4 " +3 ( 4 - 78 ) 

Nga d+1 v ' 

• (3 = 4: 

~ 2N7rga d+1 a 4n+1 (4.79) 
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5 Conclusions 



Therefore, we have obtained some exact integral expressions and asymptotics for the 
skew-orthogonal polynomials involved in the Orthogonal and Symplectic random ma- 
trix ensembles. 

Our asymptotics were derived in the "one-cut" case only, though it seems obvious 
that the result could be extended easily to the multicut case, following the method 
of [15] or |n| , it would involve hyper-elliptical theta- functions instead of exponentials. 



Another possible extension of the method presented here is to " multi-matrix mod- 
els", and a time dependant matrix, as in [ I3[ . It seems that the same kind of asymp- 
totics could be obtained. 

The asymptotics of the skew-orthogonal polynomials are useful to evaluate the 
kernels: 

N-l 

K(X, p) = V (P 2n (A)P 2n+1 (/i) - P 2n+1 (A)P 2n (/i)) e~ NV ^ e~ NV ^ (5.80) 
which give all the correlation functions. For instance with (3 = 4, we have || 



(5.81) 



p c (A,/i) = -—K(\, f x)—K(X, t i)+K(\, f ,)——K(\, t x) (5.82) 



In order to use the asymptotics of the orthogonal polynomials in p.8U| , one needs a 
generalization of the Darboux-Christoffel theorem, which allows to write K(\, fx) in 
terms of a few P n only with n close to the Fermi level np. With asymptotics of the 
type [4.66| , |4.67| , ]4.70| , (4.71j or (4.73] , the Darboux-Christoffel Theorem merely amounts 



to a formal resummation of the geometrical series (it was proven in [13] for hermitian 
multi-matrix models, and we cannot see any reason why the same proof would not work 
here). For instance in the (3 = 4 case, the generalization of the Darboux-Christoffel 
theorem reads: 

N-l 

2i sin (0(A) -M) [b * 6) 

This trick allows to find asymptotics for the kernels K(X,fi), and then asymptotics 
for all the correlation functions. One can then easilly check that in the short distance 
regime [A — y\ ~ 0(1/N), the universal 2-point connected correlation function is well 
reproduced , and that in the long distance regime |A — ~ 0(1), the smoothed 2-point 
connected correlation function is correctly reproduced too. The leading behaviour of 
short and long distance correlation functions was already known from other methods Q , 
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so that our method does not provide any new result for the correlation functions. 
However, it seems that our asymptotics can be used to build a rigorous mathematical 



proof of the universality, following the method of [10], because they allow a good control 
of the approximations. 

In addition, the fact that the skew-orthogonal polynomials are exactly the average 
characteristic polynomials of the random matrices is remarkable. It would be inter- 
esting to understand the generality of this result, and for instance try to generalize it 
to the other random matrix ensembles related to Cartan's classification of symmetric 



spaces 17, 18 
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